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ABSTRACT 



^0 ■ Context. Based on the single degenerate (SD) scenario, a super-Chandrasekhar mass model derived from the rapid rotation of a progenitor star 
' may account for the over-luminous type la supernovae (SNe la) like SN 2003fg. Chen & Li (2009) calculated a series of binary evolution and 
O showed the parameter spaces for the super-Chandrasekhar mass model. Pakmor et al. (2010) developed a equal-mass double degenerate (DD) 
I ■ model to explain sub-luminous SNe la like SN 1991bg. But they did not show the evolution of the birth rate of these peculiar SNe la or did not 
O ' compare them with absolute birth rates from observations. 

, Aims. We want to show the evolution of the birth rates of these peculiar SNe la based on their results, and compare the birth rates with 
^ ■ observations to check whether these model may account for all the peculiar SNe la. 
i 1 ' Methods. We carried out a series of binary population synthesis calculations and considered two treatment of common envelope (CE) evolution, 
i.e. a-formalism and y-algorithm. 

' Results. We found that the evolution of birth rate of these peculiar SNe la is heavily dependent on how to treat the CE evolution. The over- 
ly luminous SNe la may only occur for a-formalism with low CE ejection efficiency and the delay time of the SNe la is between 0.4 and 0.8 
Gyr. The upper limit of the contribution rate of the supernovae to all SN la is less than 0.3%. The delay time of sub-luminous SNe la from 
equal-mass DD systems is between 0.1 and 0.3 Gyr for a-formalism with a = 3.0, while longer than 9 Gyr for a = 1.0. The range of the delay 
time for y-algorithm is very wide, i.e. longer than 0.22 Gyr, even as long as 15 Gyr. The sub-luminous SNe la from equal-mass DD systems 
may only account for no more than 1% of all SNe la observed. 
CD Conclusions. The super-Chandrasekhar mass model of Chen & Li (2009) may account for a part of 2003fg-like supernovae and the equal- 
' ^ . mass DD model developed by Pakmor et al. (2010) may explain some 1991bg-like events, too. In addition, based on the comparison between 
^— ' theories and observations, including the birth rate and delay time of the 1991bg-like events, we found that the y-algorithm is more likely to 
be an appropriate prescription of the CE evolution of DD systems than the a-formalism if the equal-mass DD systems is the progenitor of 
1991bg-like SNe la. 
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1. Introduction (Zhao et a H2010l ). However, the nature of SNe la is still un- 

clear, especially its progenitor system (Hillebrandt & Niemeyer 
As one of distance indicators, type la supernovae (SNe la) Leibundgut[2000l>. It is widely believed that a SN la is 

showed their importance in determining cosmological param- produced from the thermonuclear runaway of a carbon-oxygen 
eters, which results in the discovery of_the accelerating ex- white dwarf (CQ WD) in a binary system (Arnett rT983. Based 
pansion of the universe (Riess et al. [1998J Perlmutter et al. Qn me nature of ±e WD companiolli two basic scen arios have 
[1993. The result is exciting and suggesting the presence of been discussed for about mre e decades. One is the single de- 
dark energy. At present, SNe la are proposed to be cosmologi- generate (SD) model (Whe i an & Iben [19731 Nomoto et al. 
cal probes for testing the evolution of the dark energy equation rjgg^ ; e the companion may be a main-sequence (MS) or 
of state with time andtesting the evolutionary history of the a slightly evolyed star (WD+MS) or a red giant star (WD+RG) 
universe (Riess et al. [20071 Kuznetsova et al. [200S Howell et Qr a helium ^ (WD + Hg ^ and me model h widdy in _ 
al. [20091. They even were chose to check the general relativity vestigate d by many groups (Yung elson et al. [19951 Li et al. 

[19971 Hachisu et al. ll999a|[l999bl Nomoto et al. [19991 120031 
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Langer et al. 120001 Han & Podsiadlowski 120041 Chen & Li 
120071 Han 120081 Meng & Yang 1201 Obi Lu l2009l Wang et al. 
I2009al I2009bl Wang, Li & Han |20l0l . the other is the dou- 
ble degenerate (DD) model (Iben & Tutukov 119841 Webbink 
|1984| l, in which a system of two CO WDs loses orbital angu- 
lar momentum by means of gravitational wave radiation and 
finally merges. Although the birth rate and the distribution 
of its delay time (DDT, delay time: between the episode of 
star formation producing progenitor systems and the occur- 
rence of SNe la) from this channel is comparable to obser- 
vations (Han [19981 Yungelson & Livio [19981 [2000 ; Tutukov 
& Yungelson 2002; Mennekens et al. 2010) and some SNe la 
such as SN 2003fg and SN 2005hj may likely result from the 
channel (Howell et al. 2006 1 Quimby et al. 120071) . some nu- 
merical simulations showed that the most probable fate of the 
coalescence is an accretion-induced collapse and, finally, neu- 
tron star formation (see the review by Hillebrandt & Niemeyer 
2000 ). Generally, the WD mass at the moment of explosion for 
SD model is equal (the so-called Chandrasekhar mass model) 
or less than the Chandrasekhar limit (the sub-Chandrasekhar 
mass model), while larger than the limit for DD model (super- 
Chandrasekhar mass model). 

The most remarkable property of SNe la is their apparent 
homogeneour nature. This characteristic is interpreted as the 
observable consequence of a model where the progenitors of 
SNe la are CO WDs which have increase their mass to close 
to the Chandrasekhar limit via mass transfer in a binary sys- 
tem. Nevertheless, observations of two peculiar events, i.e. SN 
1991T and SN 1991bg, raised questions about the uniformity 
of SNe la. Motivated by these discoveries, Phillips ( 1993 ) sug- 
gested a relation between the absolute magnitude of SNe la at 
maximum light and the magnitude drop of B light curve during 
the first 15 days following maximum. When SNe la are applied 
as distance indicator, the relation is adopted. 

The peak brightness of SN 1991bg was ~ 2.0 mag fainter 
in B than that of normal SNe la, and its light curve declined 
unusually fast following maximum (Filippenko et al. 1992b, 
Leibundgut et al. [T993l . SN 1991bg is identified by Branch & 
Miler (119931 ) as intrinsically subluminous. The low peak lu- 
minosity indicated a low amount of 56 Ni, i.e. about 0.07 M Q 
(Mazzali et al. I1997I ). Generally, the subluminous property of 
the SN 1991bg may be explainedby a sub-Chandrasekhar mass 
model from a CO WD + He star system or a merger of CO WD 
and He WD (Branch et al. [19951 Livio [T999ll2003l . However, 
Pakmor et al. (2010) recently suggested an alternative scenario, 
i.e. sub-luminous 1991bg-like events may be from the mergers 
of equal-mass CO WDs of M ~ O.9M . Although the light 
curve from the merger model is broader than that of 1991bg- 
like events, the synthesized spectra, red color and low expan- 
sion velocities are all close to those observed for 1991bg-like 
events. They claimed that the events from the merger of equal- 
mass WDs should occur with a rate of « 2 - 1 1% of all SNe 
la. However, Pakmor et al. ([2010) discussed a relative birth rate 
to various SNe la progenitor channels rather than the absolute 
ones from observations. One of the purposes in this paper is 
to present the evolution of the birth rate of the sub-luminous 
1991bg-like events from the mergers of equal-mass WDs and 
to compare the birth rate with those from observations. 



SN 199 IT is an overluminous event and its light curve 
declined slowly after maximum luminosity (Filippenko et al. 
I1992al Phillips et al. [19921 [19991 . which is often taken as an 
indication of a large 56 Ni mass (Hofiich et al. 1995, Nugent et 
al. 119971 Pinto & Eastman 2001). Nevertheless, Kasen et al. 
(2004) suggested a second, physically very different route to 
explain the peculiarities of SN 199 IT — one could be peering 
down an ejecta hole which is due to the existence of the com- 
panion at the moment of supernova explosion. This sugges- 
tion was uphold by detailed binary population synthesis (BPS) 
study (Meng, & Yang |2010ai) . So, it is possible that 1991T-like 
supernovae have not any special properties in physics except 
for the viewing angle of an observer. 

The SN 2003fg was observed to be 2.2 times over-luminous 
than a normal SN la and the amount of 56 Ni was inferred to 
be 1.3M G , which requires a super-Chandrasekhar mass explo- 
sion (~ 2.1M Q , Astier et al. l20"06l Howell et al. [20061 . Since 
SN 2003fg, another three 2003fg-like events were discovered, 

1. e. SN 2006gz (Hicken et al. 120071 SN 2007if (Scalzo et al. 
[20101 Yuan et al. [20T0b and SN 2009dc (Tanaka et al. l20T0l 
Yamanaka et al. 120101) . These supernovae are usually assumed 
from the mergers of DD systems, where the total mass of the 
DD systems is over the Chandrasekhar mass limit. However, a 
super-Chandrasekhar WD may also exist in a SD system, where 
the massive WD is supported by rapid rotation (Yoon & Langer 
2005; Howell et al. 120061 ). For example, Maeda & Iwamoto 
(2009) showed that the properties of SN 2003fg are more con- 
sistent with aspherical explosion of a super-Chandrasekhar WD 
supporting by a rapid rotation. Taking account of the influence 
of rotation on accreting WD, Yoon & Langer (2004) found 
that a WD under a special condition may rotate differentially. 
Using the result of Yoon & Langer d2004l Chen & Li d2009b 
calculated the evolution of close binaries consisting of a CO 
white dwarf and a normal companion, and obtained the param- 
eter space in an orbital period - secondary mass (log P;, M^) 
plane for super-Chandrasekhar SNe la. However, they also 
did not present the evolution of the birth rate of the super- 
Chandrasekhar SNe la. In this paper, we will show the evo- 
lution by a detailed BPS study. 

In Sect. [2l we describe our BPS method. We show the BPS 
results in Sect. [3] and the discussions and conclusions in Sect. 

m 

2. Binary population synthesis 

2.1. Super-Chandrasekhar SNe la 

The study of Chen & Li J2"0"09l ) is based on a WD + MS sys- 
tem and their results were summarized in an orbital period 
- secondary mass (logP^Mj) plane. According to the situa- 
tion of the primary in a primordial system at the onset of the 
first Roche lobe overflow (RLOF), there are three evolution- 
ary channels contributing to the WD + MS system, i.e. He star 
channel, EAGB channel and TPAGB channel (see Meng, Chen 
& Han [20091 and Meng & Yang 1201 Obi for details about the 
channels). In this paper, the TPAGB channel is the dominant 
one since the initial mass of WD for over-luminous SNe la is 
larger than 1.0 M , and the other two channels do not con- 
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tribute to the over-luminous because the initial WD mass from 
this channel is always smaller than 1.0 M Q (see also Fig. 8 in 
Meng, Chen & Han 2009). For all the three channels, a com- 
mon envelope (CE) phase is expected. After the formation of 
the WD + MS systems, the systems continue to evolve and the 
secondaries may also fill their Roche lobes at a stage and RLOF 
starts. We assume that if the initial orbital period, PV,, and the 
initial secondary mass, Ml, of a SD system locate in the appro- 
priate regions in the (log P\ Ml) plane for super-Chandrasekhar 
SNe la at the onset of RLOF, a super-Chandrasekhar SN la is 
then produced. 

2.2. Sub-Luminous SNe la 

Pakmor et al. ( 1201 Oi l also made a BPS study about 1991bg- 
like events based on their merging model from equal-mass DD 
systems, where the components are both massive and the mass 
ratio is only slightly less than one, i.e. the primary masses are in 
the range of 0.85 - 1 .05 M Q with a mass ratio of 0.9 < q < 1 .0. 
Since the results from the large WD mass (as large as 1 .05 M ) 
will more likely lead to more luminous SNe la, not 1991bg-like 
events, we critically limit that the range of WD mass is between 
0.83 and 0.9 M in this paper (see also Pakmor et al. 1201 Oi l. A 
primordial binary system may become a DD system after one 
or two CE phases (Han [19981 Ruiter et al. [2D09l >. There are 
three sub-categories to form the DD system contributing to the 
sub-luminous SNe la. In the following, we briefly outline the 
three sub-categories: Case 1, 2 and 3. 

Case 1 (2RLOF + 2CE): Although the evolutionary chan- 
nel for Case 1 is similar, the parameters of the primordial bi- 
naries for different CE treatment are different (see next subsec- 
tion). For a-formalism, the primordial zero-age main sequence 
(ZAMS) mass of primaries is in the range from 5.2 M Q to 5.7 
M Q , and the mass ratio {milm\) is between 0.59 and 0.7, and 
the primordial separation is close, i.e. from 30 R to 60 R . 
For y-algorithm, the primordial ZAMS mass of primaries is in 
the range from ~ 5.35 M Q to 7.1 M , and the mass ratio is 
between 0.35 and 0.53, and the primordial separation is even 
more close, i.e. from 17 R Q to 30 R Q . Due to the close sep- 
aration, the primary fills its Roche Lobe when it is crossing 
Hertzprung Gap (HG), and then a stable RLOF occurs. The 
primary loses its hydrogen-rich envelope and then becomes a 
helium star, where mass ratio reverses at a point. The helium 
star may fill its Roche lobe again after its central helium is ex- 
hausted, where the RLOF is also dynamically stable for low 
mass ratio. During the two stable RLOF, the separation may 
increase greatly for material transferring from light component 
to heavy one, and then the system finally becomes a very wide 
WD + MS system, where the separation may be as large as 400 
R Q . The secondary fills its Roche lobe during HG, which leads 
to a CE for a high mass ratio (~ 10). After the CE ejection, 
the mass donor becomes a helium star and continues to evolve. 
The helium star fills its Roche lobe once again after the exhaus- 
tion of the central helium, and a second CE ensues (mass ratio 
~ 2). Thus a DD system forms. However, the final separation 
is heavily dependent on the treatment of the CE, i.e. < Q.6R 
for the a-formalism, while 2 - 3 R Q for the y-algorithm, which 



is the dominant factor resulting in different delay time of SNe 
la for the different treatment of the CE. 

Case 2 (RLOF + CE): For the case, the mass ratio for pri- 
mordial binary system is very close to 1, i.e. larger than 0.94, 
and the primordial ZAMS mass of the primary is between 3.75 
and 4.1 M . The primordial system has a very wide separation 
(wider than 2500 R Q ), which permits the primary to evolve to 
thermal pulsing asymptotic giant branch (TPAGB) stage before 
it fills its Roche lobe, and then a stable RLOF occurs. At this 
stage, the secondary is a horizontal branch (HB) star (central 
helium burning). The system consists of a CO WD and a HB 
star after the RLOF. The secondary fills its Roche lobe again 
at TPAGB stage and a CE forms for a high mass ratio(~ 7). 
After the CE ejection, the orbit decays greatly and a DD sys- 
tem forms with a separation of 0.2 to 2 R Q . Then, the delay time 
of the sub-luminous SNe la from this sub-category is less than 
2 Gyr. 

Case 3 (one CE): The primordial ZAMS mass of primaries 
is around 4M and the mass ratio is almost 1, i.e. larger than 
0.999. The primordial separation is from 12OQR to 16OO/? , 
and thus both primary and secondary are at TPAGB stage when 
primary fills its Roche lobe, which results in a CE. After the CE 
ejection, a DD system forms with a separation of ~ 3R G . The 
delay time of SNe la from this sub-category is longer than 9 
Gyr. 

Following the DD system, gravitational wave radiation 
(GW) dominates the evolution of system on a timescale ?gw 
(Landau & Lifshitz [1962b , 

row^SxlO ^J^V , (!) 

M\M 2 

where P is the orbital period of the DD system in hours, and 
Mi and M 2 are the mass of the two white dwarf in solar mass, 
respectively. Then, the time elapsed from the birth of primor- 
dial binary system to the occurrence of SN la is equal to the 
sum of the timescale on which the secondary star becomes a 
WD and the orbital decay time. We assume that if Mi and M 2 
are both in the range of 0.83 - 0.9 M and the elapsed time is 
less than 15 Gyr, a sub-luminous SN la is produced. 

2.3. Common envelope 

As mentioned above, CE is very important for the formation 
of SD and DD systems. During binary evolution, the mass ra- 
tio (q = Mdonor/^accretor) is crucial. If it is larger than a critical 
mass ratio, q c , mass transfer between the two components is dy- 
namically unstable and a CE forms (Paczyriski 1976). The ratio 
q c varies with the evolutionary state of the donor star at the on- 
set of RLOF (Hjellming & Webbint fi987l Webbink [l988l Han 
et al. 120021 Podsiadlowski et al. l2002l Chen & Han [2008]). In 
this study, we adopt q c = 4.0 when the donor star is on MS or 
crossing HG. This value is supported by detailed binary evolu- 
tion studies (Han et al. l2000l Chen & Han l2002ll2003l ). If the 
primordial primary is on first giant branch (FGB) or asymptotic 
giant branch (AGB), we use 

q c = [1.67-x + 2(^) 5 ]/2.13, (2) 
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where M c is the core mass of the donor star, and x = 
dln/?i/dlnM is the mass-radius exponent of the donor star 
and varies with composition. If the mass donors (primaries) are 
naked helium giants, q c = 0.748 based on Equation (O (see 
Hurley et al. l20"02l for details). 

Embedded in the CE are the dense core of the donor star 
and the secondary. Due to frictional drag with the envelope, the 
orbit of the embedded binary decays, and a large part of the 
orbital energy released in the spiral-in process is injected into 
the envelope (Livio & Soker 1988). It is generally assumed that 
the CE is ejected if 



(3) 



where AE m \, is the orbital energy released, £bind is the binding 
energy of CE, and uqe is CE ejection efficiency, i.e. the frac- 
tion of the released orbital energy used to eject the CE. Since 
the thermal energy in the envelope is not incorporated into the 
binding energy, q-ce may be greater than 1 (see Han et al. 1995 
for details about the thermal energy). In this paper, we set acE 
to 1.0 or 3.CQ. After the CE, the orbital separation at the onset 
of the CE, a;, becomes Of which is determined by 



G(M C + M e )M e 
AR, 



GM c m GMm 
«ce( — ), 



2flf 



(4) 



where A is a structure parameter relaying on the evolutionary 
stage of the donor, M, M c and M e are the masses of the donor, 
the donor envelope and the core, respectively, Ri is the radius of 
the donor, and m is the companion mass. In the paper, we take 
the structure parameter a constant, i.e. A = 0.5 (de Kool, van 
den Heuvel & Pylyser 1987). Thus, the final orbital separation 
of a binary system after the CE phase cif is given by 



Of _ Mc 



2M e ai 
acE^tnR 



■r 1 . 



(5) 



The description above may well produce the distribution of 
WD + MS systems as noticed by Zorotovic et al. (2010) (see 
also Hurley et al. 2002), while it is possibly difficult for the de- 
scription to produce a close pair of white dwarfs. Nelemans et 
al. ((2000) and Nelemans & Tout ( 2005 ) suggested an alterna- 
tive algorithm equating the angular momentum balance, which 
may explain the formation of all kinds of close binaries: 



A/ 
T 



tce 



M + m 



(6) 



where J is the total angular momentum and AJ is the change 
of the total angular momentum during the common envelope 
phase. Implicitly assuming the conservation of energy, the final 
orbital separation Of after the CE is then given by 



Of M 2 M c + m 
^ ~ Mc M + m 



)(1 - Tce 



Me 



M + m 



(7) 



Based on the results in Nelemans & Tout (2005), we set jqe 
1.5. 



1 Via reconstructing the evolution of post-common-envelope bina- 
ries (PCEB) consisting of a white dwarf and a main-sequence star, 
Zorotovic et al. (2010) obtained a relative low oce value, where the 
internal energy of the envelope was included. 
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Fig. 1. Evolution of the birthrates of super-luminous SNe la 
from SD channel and sub-luminous SNe la from the merge of 
equal-mass DD systems for a single starburst of 10 u M o for 
different CE treatments. The thick line is the result for over- 
luminous SNe la, while the thin lines are the results for sub- 
luminous SNe la 



Following Nelemans & Tout (2005), we call the the formal- 
ism for equation Q the a-formalism and that of equation ((6} 
the y-algorithm. In this paper, we treat CE evolution by the two 
algorithms. 

2.4. Basic parameters for Monte Carlo simulations 

To investigate the birthrate of SNe la, we followed the evolu- 
tion of 4 x 10 7 binaries via Hurley's rapid binary evolution code 
(Hurley et al. 2000, 2002). The descriptions above for super- 
luminous SNe la from the SD model and the sub-luminous 
ones from the equal-mass DD systems are incorporated into the 
code. The primordial binary samples are generated in a Monte 
Carlo way and a circular orbit is assumed for all binaries. The 
basic parameters for the simulations are as follows: (1) a con- 
stant star formation rate (SFR) of 5 M yr _1 over the past 15 
Gyr or a single star burst of 10 n M G ; (2) the initial mass func- 
tion (IMF) of Miller & Scalo ( (19791 : (3) the mass-ratio dis- 
tribution is constant; (4)all stars are members of binary sys- 
tems and the distribution of separations is constant in log a for 
wide binaries, where a is the orbital separation, and falls off 
smoothly at small separation, where a = 10R Q is the boundary 
for wide and close binaries; (5) solar metallicity, i.e. Z = 0.02 
(see Meng & Yang (2010b ) for details of the parameter input). 

3. Results 

3.1. Single star burst 

Fig. Q] shows the evolution of the birth rate of over-luminous 
(the thick line) and sub-luminous (the thin lines) SNe la with 
different CE treatments, where a single star burst is assumed. 
An imposing property in the figure is that the results are heavily 
dependent on the treatment of CE evolution. For over-luminous 
SNe la, only a-formalism with or = 1 .0 may produce SNe la. 
This is mainly derived from constraint of the parameter space 
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Table 1. The Galactic birth rate for super- (Row 3) and sub-luminous (Row 4) SNe la with different BPS parameters (Row 2), 
where a constant star formation rate (5 M yr _1 ) is assumed over the past 15 Gyr. 



Series Number 




parameters 




v^lO^yr 1 ) 


v 2 (10- 6 yi'~') 


1 


«CE 


= 1.0; n(q) 


= 1; e = 





7.98 


0.359 


2 


<*CE 


= 3.0; n(q) 


= l;e = 





0.0 


12.7 


3 


7CE 


= 1.5; n(q) 


= 1; e = 





0.0 


32.6 
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ffCE 


= 1.0; n(q) 


= 2q; e - 
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1; n{e) 


= 2e 
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Fig. 2. Similar to Fig.Q]but for a constant star formation rate of 
SMoyr 1 . 

for over-luminous SNe la and the different evolutionary chan- 
nels for WD + MS systems. Based on the results in Chen & 
Li ([2009 ), a CO WD may explode as an over-luminous SN la 
only when its initial mass is larger than 1.0 M , As described 
in section |2~T1 there are three channel may produce WD + MS 
system, but the TPAGB channel is the only one leading to the 
WD + MS system with an initial WD mass larger than 1.0 M . 
The primordial binary system experiencing the TPAGB chan- 
nel has a very long orbital period, i.e. longer than 1000 days, 
which means a low binding energy of CE (Meng, Chen & Han 
2009). Because of the low binding energy of the CE and a 
long primordial orbital period, acE has a remarkable influence 
on CO+WD systems from the TPAGB channel. Generally, if 
a CE can be ejected, a low acE produces a shorter orbital- 
period WD + MS system, which is more likely to fulfill the 
conditions for over-luminous SNe la. Therefore, we only ob- 
tain over-luminous SNe la when acE = 1-0 (see also Meng, 
Chen & Han 2009). For WD + MS systems, the effect of the y- 
algorithm is similar to that of the a-formalism with acE = 3.0 
(Nelemans & Tout 2005). Thus no over-luminous SNe la is 
produced for the y-algorithm. 

It is clearly shown in Fig.Q]that most of over-luminous SNe 
la occur between 0.4 and 0.8 Gyr. This is mainly determined 
by the high initial companion mass of WD (see Figs. 2 and 3 
in Chen & Li 2009). In addition, the peak value of the birth 
rate for over-luminous SNe la is much smaller than that for 



normal SNe la from WD + MS systems by about two orders 
of magnitude (see Han & Podsiadlowski 2004; Meng, Chen & 
Han |2009l . 



Both tt-formalism and y-algorithm may produce the sub- 
luminous SNe la. However, the sub-luminous SNe la may be 
from different sub-categories for different CE treatment. For 
ff-formalism, sub-luminous SNe la are produced from Case 1 
only when acE has a high value, i.e. acE = 3.0. Because a bi- 
nary system have experienced two CE phase for Case 1 before 
it become a DD system suitable for sub-luminous SNe la, a 
low acE means a shorter orbital period, even a merger before 
CE ejection, and then no sub-luminous SNe la. Actually, when 
acE = 1 -0, all the binary system suitable for sub-luminous SNe 
la merger before a DD system forms, except those from Case 
3. Since the primordial binary systems from Case 3 have a very 
large separation, its final separation is also very large after a 
CE phase. Thus, the delay time from this channel is very long, 
i.e. longer than 9 Gyr. Whatever, the birth rate from this chan- 
nel is too low, i.e. the peak value is only 10~ 6 yr~' (see the 
thin solid line in Fig.Q}. Even for acE = 3.0, the result of sub- 
luminous SNe la significantly differ from that for y-algorithm. 
The sub-luminous SNe la occur between 0.1 and 0.3 Gyr for 
a = 3.0, while explode with a very wide age range for y- 
algorithm, i.e. longer than 0.22 Gyr and even as long as 15 
Gyr. The sub-luminous SNe la for y-algorithm come from two 
sub-categories, i.e. Case 2 with delay time shorter than 2 Gyr, 
while Case 1 with delay time longer than 2 Gyr. 



Another imposing property in Fig.Q~]is that the peak value 
of the birth rate is much lower than those from observations by 
at least two orders of magnitude (Totani et al. 2008 ; Maoz et al. 
I2010al Maoz & Badenes l20T0t . and the shape of the DDT in 
this figure also significantly differs from observations (Maoz et 
al. l2010bV The best fitted DDT from observations follows t 
(Totani et al. 120081 Maoz et al. 1201 Obi ), while the DDT is a nar- 
row peak for the super-luminous SNe la from the SD channel 
and the sub-luminous ones from the merge of equal-mass DD 
systems under <?-formalism. Although the DDT for the sub- 
luminous under the y-algorithm is not a narrow peak, its shape 
still does not follow t~ x . These results are good evidence for 
sub-group nature of 2003fg-like and 1991bg-like SNe la. 
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3.2. Constant star formation rate 

Fig. [2] shows Galactic birthrates of SNe la for over-luminous 
and sub-luminous SNe la. Our simulations give an upper limit 
of the birth rate of over-luminous SNe la from WD + MS sys- 
tems, i.e. ~ 8 x 10~ 6 yr _1 (thick solid line), which is smaller 
than the Galactic birthrate inferred observationally by about 
three orders of magnitude (3-4xl0~ 3 yr~', van den Bergh & 
Tammann [19911 Cappellaro & Turatto [19971 Li et al. l20T0l l. 
So, the over-luminous SNe la from WD + MS systems could 
be very rare events. The upper limit of the birth rate of sub- 
luminous SNe la from the equal-mass DD systems is only 
slightly higher than that of over-luminous SNe la, i.e. ~ 3 x 
10~ 5 yr~', and then the sub-luminous SNe la from DD systems 
are also rare events. 

4. Discussions and conclusions 

4.1. Birth rate 

According to the results of Chen & Li ([2009), we found that 
the 2003fg-like supernovae from SD systems contribute to no 
more than 0.3% of all SNe la. The contribution rate is so low 
and seems not to account for the discovery of four 2003fg-like 
events at present. Considering the delay time of the 2003fg-like 
supernova from the SD systems should be less than 0.8 Gyr 
while SN 2009dc may be from an old population, the model of 
Chen & Li (2009 ) may thus account for a part of the 2003fg- 
like supernovae. 

In this paper, we found that the sub-luminous SNe la 
from equal-mass DD system may at most account for 1 % 
of all SNe la observed. But Pakmor et al. (120101 ) found that 
their model should occur with a rate of 2 - 11% of the 
total SN la rate, which seems consistent with the observation 
rate of 1991bg-like supernovae with errors (16% + 7%, Li et 
al. 120011 . However, their high birth rate is relative to other 
possible SNe la progenitor formation channel, rather to that 
from observations. Since no progenitor model may account for 
the birth rate of SNe la derived observationally at present, 
the birth rate in Pakmor et al. (120101 ) might be overestimated. 
Then, the equal-mass DD systems may be the progenitors of 
somel991bg-like supernovae. 



4.2. Age 

In the paper, we find that if an over-luminous SNe la is from 
a WD + MS system, its delay time should be shorter than 0.8 
Gyr, which means that there is new star formation during the 
recent 0.8 Gyr in the host galaxies of the over-luminous SNe 
la if the SNe la are from WD + MS systems. At present, 
there are four documented examples of over-luminous SNe 

2 The DD model may explain the observational birth rate of SNe 
la (Han 1998; Yungelson & Livio 1998 2000), but earlier numerical 
simulations showed that the most probable fate of the coalescence of 
DD systems is an accretion-induced collapse and, finally, neutron star 
formation (Saio & Nomoto 1985, see also the review by Hillebrandt 
& Niemever l2000l 



Ia explosions with an derived progenitor mass higher than 
Chandrasekhar mass limit, i.e. SN 2003fg, 2006gz, 2007if and 
2009dc. Amoung the four SNe Ia, only the host galaxy of SN 
2009dc was suspected to have no significant information of star 
formation because it is an SO galaxy and inspection of its SDSS 
images show uniformly red color (Yamanaka et al. 120 101) . 

The delay time of 1991bg-like events is much different 
based on the treatment of CE. For a = 3.0, the delay time is 
shorter than 0.3 Gyr, while longer than 9 Gyr for a - 1.0. 
The range of delay time for y-algorithm is very wide, longer 
than 0.22 Gyr and even as long as 15 Gyr. Thus, our results 
provides a possibility to test which treatment of CE may work 
for the CE evolution of equal-mass DD systems via measuring 
the progenitor age of 1991bg-like events. Now, many SNe Ia 
like SN 1991bg were discovered and it is apparent that they 
favor E and SO galaxies, i.e. more than 60% of 1991bg-like 
events are found in E/S0 galaxies and the remainder are found 
in early-type spiral galaxies (Howell 2001 ), and then the sub- 
luminous SNe Ia are expected to arise in old (> 1 Gyr) stellar 
population (Sullivan et al. 2006). For example, the host galaxy 
of SN 1991bg, NGC 4371, is an elliptical galaxy and its age 
is approximately 14 + 2.5 Gyr (Trager et al. 120001) . Although 
this age relies on the assumption of a single star burst and is 
model dependent, the results of Trager et al. ((2000) provided 
a lower limit of 10 Gyr for the age of the progenitor of SN 
1991bg (Howell 2001). The simulation of a - 1.0 seems to 
favor such long age, but the birth rate from the simulation is 
too low to compare with observations, and the age range for 
the simulation is also too narrow. Considering that the simula- 
tion of a = 3.0 also can not explain the old population nature 
of 1991bg-like SNe Ia, the a-formalism may not explain the 
properties of these sub-luminous events. In our simulation for 
y-algorithm, the Case 1 (older than 2 Gyr) dominates the sub- 
luminous SNe Ia and contributes to more than 95% of these 
SNe Ia. This means that most of 1991bg-like events belong to 
old population, while there is still a few the sub-luminous SNe 
Ia favoring young population, which is more consistent with 
the statistic results of Howell (2001) than a-formalism, at least 
qualitatively. In addition, the birth rate for the y-algorithm is 
also more close to that from observations (Sullivan et al. 2006). 
So, our results are more likely to favors the y-algorithm[| One 
may argue that since the simulation of a = 1 .0 may produce the 
old sub-luminous SNe Ia as obtained in Pakmor et al. (1201 Oi l, 
my conclusion is possibly dependent on the parameters of BPS 
simulation, such as the distributions of mass ratio and eccen- 
tricity. We check the influence of the distributions of mass ratio 
and eccentricity on the birth rates of the special SNe Ia under 
the assumption of a-formalism. We find that the influence on 
the super-luminous SNe Ia is insignificant (see TableQ]), while 
on the sub-luminous SNe Ia can not be neglected. The distri- 
bution of n{q) = 2q favors the high mass ratio, and then may 
enhance the contribution of Case 3 and the birth rate of the sub- 



3 Note: after analyzing the evolution of a new sample of PCEBs, 
Zorotovic et al. (2010) found that the classical a-formalism seems 
to be an appropriate prescription of CE evolution and turns out to 
constrain the outcome of the CE evolution much more than the y- 
algorithm. 
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luminous SNe la, but the enhancement is moderate by about a 
factor of 2 and the birth rate is still much lower than that derived 
observationally (see Table [TJ. Elliptical orbit is more possible 
to deduce the CE evolution than circular orbit. So, under the 
distribution of n(e) = 2e, the contribution of Case 1 is signif- 
icantly enhanced, and then the birth rate of the sub-luminous 
SNe la (see TableQ]i. However, the distribution of n(e) = 2e al- 
most dose not affect Case 3. Since the delay time from Case 1 is 
very short, i.e. shorter than 0.3 Gyr, most of 1991bg-like events 
(more than 97%) should be from very young population if a- 
formalism works for the sub-luminous SNe la, which conflicts 
with observations (Howell 2001). So, our basic conclusion can 
not be influenced by BPS parameters^. However, please keep in 
mind that our conclusion is critically dependent on the result of 
Pakmor et al. (2010), i.e. the 1991bg-like supernovae are from 
equal-mass DD systems, which are only a part of SNe la from 
DD systems, not all DD systems. 

In summary, we show the evolution of the birth rate of over- 
luminous SNe la from SD systems and sub-luminous SNe la 
from equal-mass DD systems in this paper, and found that the 
birth rate of the SNe la are very low. Our results indicate that 
the model of Chen & Li (2009 ) may account for some 2003fg- 
like supernova, so as to the model of Pakmor et al. (20103 for 
1991bg-like events. In addition, depending on the simulation 
of Pakmor et al. (|2010), the results in the paper favors the y- 
algorithm for the CE evolution of equal-mass DD systems. 
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